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ABSTRACT 


An attempt is made to estimate aerodynamic parameters using radar tracked 
coordinates of a free flight artillery shell. The present thesis proposes schemes to estimate 
average values of aerodynamic parameters and their functional dependence on Mach 
number using Maximum Likelihood method with limited number of instrumented data. 
The suitability and applicability of the proposed schemes are evaluated by applying them 
to 2-dimensional and 3-dimensional flight data. 

The proposed schemes were firstly validated using 2-dimensional simulated flight 
data and as the work progressed, the trajectory model was made 3 -dimensional. The 
results obtained using these schemes show good accuracy of estimated parameters. The 
Cramer-Rao bounds were low even when moderate noise was included. However 
limitations of few schemes to deal with limited flight data for parameter estimation 
purpose is also highlighted. 
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NOMENCLATURE 


C L , C d0 , C m = Non-dimensional lift, drag and pitching moment coefficients. 

C N ,C A = Coefficient of normal and axial force. 


C x C y , C, = Coefficient of longitudinal, lateral, vertical force. 
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u, v,w 
V 

W x ,W y ,W z 


= Diameter of the artillery shell, mm. 

= Acceleration due to gravity, m/s 2 . 

= Moment of inertia about x and y axis, kgm 2 . 

= Mass of the shell, kg. 

= Roll rate, rad/s. 

= Ballistic Air Pressure (ambient Air Pressure), mm Hg 
= Pitch rate, rad/s. 

= Dynamic pressure, N/m 2 . 

= Yaw rate, rad/s. 

= Distance from earth center to C.G. of the shell, m. 

= Radius of earth, m. 

= Reference area of the artillery shell, m 2 . 

= Time of flight, sec. 

= Ballistic air temperature (ambient air temperature), °C. 

= Velocity components in x, y and z body axes, m/s. 

= Air relative speed, m/s. 

= Head/tailwind, crosswind and vertical wind components, m/s. 


x,y,z 


= Spatial coordinates, m. 


a r ~ Yaw of repose. 

c o = Rotation vectors. 

Q = Earth rotation angular velocity, rad/s. 


IX 



p 

0 

</> 


= Density of air, kg/m 3 . 

= Firing Table Elevation, Pitch attitude, deg. 

= Roll attitude, deg. 


V 


a 

fi 


= Yaw angle, deg. 

= Angle of attack, deg. 
= Angle of sideslip, rad. 


Superscripts 

. = Derivative with respect to time. 

-» = Vector quantity. 


Subscripts 

o 

x,y,z 

wind 


= Initial conditions. 

= Components along x, y and z direction. 
= Wind axes. 


Stability and control derivatives 

r _dC L dC L J5Q 

da ’ Lq d(qd/2Vy Ls dS ' 
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r -?£l r dCl 

15 ~ dS ’ ,p d(pd/2V)' 

r J£il r dC " 

np dp’ nr d(rd/2V)' 

dC y dC y 

Cyp ~ dp’ Cyr ~ d(rd/2V)' 



CHAPTER 1 

INTRODUCTION 


Artillery forms an important wing of the army to provide firepower during war as 
well as during cross-border skirmishes with the enemy. Artillery generally falls into three 
basic categories: guns, howitzers and mortars. The principle difference between them 
being the trajectory of the round fired. A gun has a high muzzle velocity and a very flat 
trajectory. Normally a gun is used in a direct fire mode where the target can be seen and 
penetration is desirable. Howitzers have a somewhat lower muzzle velocity and arc their 
shells onto a target. They are used in both a direct fire and indirect fire mode. This is 
especially useful when an enemy is concealed behind a prepared position or the artillery 
men desire to have a shell explode over an enemy’s head. The airburst does less damage 
to hardened targets, but causes many more human casualties due to shell fragmentation 
covering a large area. Mortars have a very pronounced arc of flight. They have a 
relatively low muzzle velocity and are unsuitable for direct fire. Their principle value 
comes from being able to lob shells behind an obstacle, such as a fortification or a hill. 
They are not very accurate and dependent upon the amount of propelling powder to 
determine the point of impact. 

The effectiveness of artillery is largely judged by the accuracy in hitting the 
targets. The accuracy and reliability is influenced by the design criteria used in designing 
the shell. Artillery shells are a class of projectiles around which much of the aeroballistic 
theory was originally developed, and they continue to form a significant part of the 
aeroballistician's interest. 
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Study of the motion of projectiles through an external medium is known by the 
name of external ballistics. If the external medium is the earth's atmosphere, then external 
ballistics become synonymous with aeroballistics. 

The conventional approach hitherto for understanding the in-flight behaviour of 
projectiles was to develop mathematical models that could predict all elements of the 
trajectory from launch to target. To this purpose, it becomes essential that all forces, 
moments affecting the flight of the projectile be accounted for in a well-defined 
mathematical form 1 . Beginning with the most simple but relatively inaccurate 
mathematical model, the in-vacuo trajectory model, more and more sophisticated models 
of increasing accuracy such as the point mass model (PMM), the modified point mass 
model (MPMM) and the six-degrees-of-ffeedom (SDF) model have been developed 1 . 
These models require aerodynamic coefficients as input (e.g., drag coefficient C d , 

damping in roll derivative C lp , lift curve slope C La , etc.). Accurate values of these 

parameters are required in mathematical model for description of artillery shell dynamics. 
There are three distinct approaches for estimating aerodynamic parameters: 

□ Theoretical methods. 

□ Wind tunnel testing. 

□ Flight testing. 

At the preliminary design stage of any system, theoretical methods 2, 3 ’ 4 are useful 
inspite of their limited accuracy. The wind tunnel testing improves the accuracy of 
estimation but it is a time consuming and expensive way of estimating the aerodynamic 
parameters. Precise simulation of control surfaces, power effects and flight conditions is 
difficult. The model tested in the wind tunnel is generally slightly different from actual 
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flight due tQ last minute configuration changes. Other reasons for discrepancies between 
flight and wind tunnel results are Reynolds number discrepancies and interferences due to 
support system. It is therefore desirable that the wind tunnel estimates be corroborated 
with the estimates from actual flight test data. 

The most commonly used methods for parameter estimation from flight data are 
broadly classified into the following categories: 

□ Equation error methods. 

□ Output error methods. 

The principle of Least squares is used in the equation error method wherein the 
error gets minimized with respect to the unknown parameters in each of the equations. Its 
advantages include computational simplicity, non-iterative nature and applicability to 
both linear and non-linear models. The disadvantage is that the method cannot be directly 
applied if all the states are not measured accurately, and produces poor results if the 
measurements are noisy. Therefore in order to get accurate results, considerable effort is 
required for data reconstruction and smoothing. Sometimes these data processing tasks 
are more complicated than the task of parameter estimation itself. However these 
methods are useful as startup for more advanced estimators. 

In the output error method, the error between the measured and the model response 
is minimiz ed. The method assumes that there exist no modeling errors. This method 
processes the measurement noise while assuming the model representation of the given 
system to be exact. A comprehensive survey of these methods is reported by Maine and 
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IliS 5 Th i 

met ods like analog matching 6 , Newton Raphson method 7 , modified Newton 
Raphson method 7 etc. fall in this categoiy. 

Advanced methods are capable of estimating parameters in the presence of 

nt and/or process noise. Maximum Likelihood (ML) estimates are those for 

which the observed value would be the most likely to occur; “most likely” is defined as 

tnaxamzation of Likelihood function of the observed responses, given the parameters. 

The 

mam advantage of this method is that parameter estimates are asymptotically 
unbiased, consistent and efficient provided the model assumed is correct. 

imation activities in aerospace science started in late fifties. Though plenty of 
literatures on parameter estimation exist for aircraft very few exist for artillery shells. It is 
mentioning that parameter estimation technique is a blend of both science and art. 
arameter estimation of an aircraft with many measurements like, acceleration (both 
linear and angular), angular orientation, speed, angle of attack etc. has been discussed in 
fr°m cost effectiveness point of view (and also because of non- 
availability of space for instrumentation), it may not be feasible to use many sensors for 
illery shells going through many development flight trials. 

As an alternative approach, limited flight tests are conducted with scale down 
odel in Aero-ballistic range facility 9 . The test facility is an enclosed, atmospheric, 
instrumented concrete structure used to examine the exterior ballistics of various free 
flight configurations. The facility contains gunroom, control room, model measurement 
room, blast chamber and instrumented range. The range may have large number of 
locations available as instrumentation sites. Each station is used to house folly 
instrumented orthogonal shadowgraph stations. Using this instrumentation, spatial 


4 



positions, speed and angular orientations etc. are estimated. These are then used for 
parameter estimation. 

Winchebach et al 9 conducted free flight aeroballistic tests to obtain subsonic and 
supersonic aerodynamic parameters of a wrap around fin configuration. The aerodynamic 
parameters presented were extracted from the position-attitude time histories of the 
experimentally measured trajectories using non-linear numerical interpolation data 
reduction routines. 

Defense Research Establishment, Quebec, Canada is currently engaged in 
research on high l/d energy penetrators for 70mm air launched rocket application. One of 
the objectives of the program is to enlarge the free-flight aerodynamic database of high 
l/d finned projectiles 10 . Dupis investigated she model configurations (scale down) during 
the test program 10 . Standard SDF routine with ML method was used to match the 
theoretical trajectory to the experimentally measured trajectory. In that analysis, as a first 
step, the basic dynamic range data (time, position, angles), physical properties and 
atmosphere conditions were assembled and then SDF (with ML) were performed to 
estimate the aerodynamic parameters 10 . The results of this analysis showed that once time 
histories of position, angles are recorded it is possible to estimate aerodynamic 
parameters using free flight test data. 

It is worth mentioning that aerodynamic parameters are strong functions of Mach 
number. In order to capture this dependence on Mach number, projectiles are fired at 
different initial Mach number and the basic data corresponding to that Mach number are 
recorded and analyzed for parameter estimation 9 ’ 10 . Though this approach looks attractive 
(if such facility is available), however it has following two distinct limitations: 
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1) It may not be possible to simulate actual flight path that the full-scale model would be 
traversing. 

2) It requires large number of firings to capture Mach number dependence of 
aerodynamic parameters. 

In view of the above, there have been conscious efforts towards generation of flight 
data of free flight projectile by simultaneous measurements with pulse radar or optical 
system for position data and continuous wave Doppler radar for velocity measurements. 
In order to improve the overall accuracy several systems of each type have often been 
used in parallel 11 . However, this complicated and very expensive set up is not required 
today. From a single continuum wave Monopulse Azimuth and Elevation Doppler 
Tracking Radar it is possible to extract highly accurate velocity and position data in three 
dimensions without any assumptions of ballistic behaviour 11 . 

Defense Research Organization of India, is currently engaged in design and 
development of new artillery shell, rockets etc. There is a routine demand for post flight 
analysis of flight data. In particular, there is strong demand for estimation of aerodynamic 
parameters from flight data to validate the design as well as to prepare database for fire 
control system (FCS) to enhance artillery effectiveness. Accurate values of aerodynamic 
parameters are of paramount importance in building effective FCS. 

Due to non-availability of Aero-ballistic range facility in India, lots of efforts are 
being directed towards acquiring flight data through Radar tracking. At present, the radar 
system available can only measure spatial distance (x ,y,z) accurately. The available 
radar does not have facility to measure velocity, especially for long range artillery shell. 
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Thus, the motivation for this research is inspired by this challenge of estimating 
aerodynamic parameters using radar tracked spatial coordinates of artillery shell. 

Free flight projectiles experience aerodynamic forces and moments as it traverse 
from launch point to the target end. Bryan 12 introduced the concept of flight vehicle 
aerodynamic modeling, which provided relationship between the forces along the three 
Cartesian axes and the three moments about these axes as a function of linear 
translational motion variables and rotational rates. The most commonly employed 
aerodynamic models are typically formulated as truncated Taylor-series expansion of the 
form. 


c -?f^ 


i = x,y, z,l,m,n. 


where X are the vehicle motion variables (angle of attack, angle of sideslip, angular 
rates etc.). 

Conventionally, application of ML method will require measurements of motion 


variables for the estimation of the aerodynamic derivates (parameters, 



However, 


for the problem chosen in this thesis, we do not have the information of these motion 
variables. The proposed method uses only the information about the spatial location of 
the shell to estimate aerodynamic parameters. 

Artillery shells are fired at high muzzle velocity (2-3 Mach) and are expected to 
reach maximum range at an angle of elevation of around 45 degrees. During this motion 
it encompasses high supersonic to low subsonic Mach numbers. The aerodynamic 
parameters are strong function of Mach number. 
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In order to develop accurate mathematical model, it is necessary to capture the Mach 
number dependence of these parameters. In order to capture their dependence, several 
schemes have been proposed to estimate aerodynamic parameters from tracked data 
(x,y,z) of an artillery shell. Through the proposed schemes, attempt has been made to 
estimate equivalent (constant) values of aerodynamic parameters that are expected to 
account for the dependence of parameters on Mach number. Further schemes have been 
proposed to capture the Mach number dependence of these parameters by modeling the 
functional form of their variation in estimation algorithm. 

Due to non-availability of the real flight data, all the schemes as stated above were 
tested on simulated flight data. To start with, simulated data was generated using PMM. 
This is a simplified model and only uses drag coefficient (C d ) for trajectory computation. 
The proposed schemes were applied to estimate equivalent as well as Mach number 
dependent of C d successfully. Next the simulated data was updated by incorporating 

three more aerodynamic parameters by generating trajectory data using MPM model. 
Proposed schemes were applied to estimate the aerodynamic parameters. Finally, the 
simulated flight data was generated using SDF model, and proposed schemes were 
applied to investigate its capability to estimate large number of aerodynamic parameters 
used in simulation. The difficulties encountered are also presented as the work concluded 
with future suggestions. 

The work carried out in this thesis is presented in 4 chapters. Chapter 1 consists of 
brief introduction about parameter estimation and Chapter 2 describes various schemes 
proposed and generation of simulated flight data for estimation. Chapter 3 consists of 
results and discussions followed by conclusions and future suggestions in Chapter 4. 
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CHAPTER 2 


PROPOSED SCHEMES AND GENERATION 

OF FLIGHT DATA 


2.1 General 

This chapter describes the various schemes proposed and different types of 
simulated flight data generated for validating the proposed schemes. Various 
mathematical models consisting of the equations of motion that had been used for the 
generation of flight data are being outlined. 

2.2 Proposed schemes 

The following schemes have been proposed for estimation of aerodynamic 
parameters from flight data of an artillery shell using ML method: 

Scheme 1: In this scheme motivation was to estimate equivalent parameters, 
which can be assumed to be an approximation for the variation of parameters with Mach 
number. In order to validate this scheme flight data were generated, where the variation 
of aerodynamic parameters with Mach number were used in the simulation model. 
However, the Mach number dependence was not modeled in estimation algorithm, and 
thus estimated parameters (equivalent parameters) were expected to include the effects of 
Mach number variation in some approximate manner. 

Scheme 2: In this scheme motivation was to estimate the functional dependence 
of parameters on Mach number. In order to validate this scheme flight data were 
generated by incorporating variation of aerodynamic parameters with Mach number (as 
supplied by the manufacturer) in simulation model. However, in the estimation algorithm 
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the Mach number dependence of aerodynamic parameters were approximated by 
expressing aerodynamic parameters as a polynomial function of Mach number. The 
constants of this polynomial series were estimated applying ML method. 

Scheme 3: In this scheme motivation was same as that of scheme 2 but different 
methodology was used. The whole trajectory was split into different sets of points and 
average values of the parameters pertaining to each set of points were estimated. The 
estimated parameters correspond to average Mach number of the different set of points 
used for estimation. 

Scheme 4: Finally, in this scheme, the flight data was generated at a particular 
elevation corresponding to different charges (muzzle velocities), and small initial part of 
the trajectory was used to estimate the average values of aerodynamic parameters. These 
estimated parameters correspond to the Mach number corresponding to different muzzle 
velocities. 

As the simulation was made exhaustive by incorporating large number of 
parameters, few of the schemes faced difficulties to estimate all the parameters. Flight 
data were generated using all the three trajectory models as explained in next section. The 
aerodynamic and geometric characteristics of the artlillery shell used for generating the 
simulated flight data are listed in Table 2.1 and Table 2.2 respectively. 
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Table 2.1 Aerodynamic characteristics of the artilleiy shell 


SL. No. 

Mach 

Number 

Cd 

c 

ma 

Q* 

C/, 

1. 


0.13 

3.2 

1.3 

WEMm 

2. 

0.85 

0.13 

3.2 

1.3 

MKjmmm 

3. 

0.9 

0.13 

3.25 

1.2 

-0.03 

4. 

ME^Mi 

0.14 

3.3 

1.1 

BEESBI 

5. 

HKE9H 


3.45 

6.95 

■EiXSOfll 

6. 


0.269 

3.5 

0.95 


7. 

i 

0.329 

3.5 

1 

— 

8. 

1.02 

0.351 

3.4 

1.1 


9. 

1.05 

0.367 

3.25 

1.55 

-0.029 


1.1 

0.366 

3.15 

1.85 

-0.029 

11. 

1.15 

0.353 

3.15 

1.95 

-0.0285 

12. 

1:2 

0.343 

3.14 

2 

-0.0282 

13. 

1.3 

0.326 

3.13 


BEEn 

14. 

1.4 


3.12 

2.15 

1 

15. 

1.6 

0.293 

2.75 

2.2 

-0.026 

16. 

1.8 

0.273 

2.5 

2.3 

-0.025 

17. 

2 

0.257 

2.48 . 

2.32 

aEESH 

18. 

2.2 


2.47 

2.35 


19. 

2.4 

Ha 

2.45 

2.35 

aEsm 


2.6 

0.229 

2.4 


-0.0215 


Table 2.2 Geometric characteristics of the artillery shell 


SL. No. 

Geometric characteristics 

1 

Total length 

: 825 mm. 

2 

CG location 

: 533 mm from nose tip. 

3 

Diameter 

: 155 mm. 

4 

Mass 

: 42.6 kg. 

nr* 

5. 

Moment of Inertia 


Axial, 

: 0.146 kg-m . 


Traverse,/^ :1.709kg-m 2 . 
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2.3 Trajectory models used for simulating flight data 

In general, the trajectory followed by center of mass of artillery shell in space is 
3-Dimensional (x,y,z). Theoretical modeling of such trajectory demands aerodynamic 
modeling to be exhaustive. It requires usage of large number of aerodynamic parameters. 
Hence estimation of all parameters at one time using proposed schemes might be 
involved. It was thus decided firstly to validate/investigate the proposed schemes using 
theoretically simulated 2-Dimensional flight data. As the work progressed, the trajectory 
model was made 3 -Dimensional and all the schemes were tried. Accordingly the work 
started using simulated data through Point Mass Model (PMM) and then extended to 
Modified Point Mass Model (MPMM) and Six Degree of Freedom model (SDF) *. 


2.3.1 Point Mass Model (PMM) 

In this model, it is assumed that the only aerodynamic force acting on the 

projectile is drag. It provides fairly accurate estimates of range for adequately stable 

projectiles and can also be used to estimate the first order effects of wind. The equations 

of motion for point mass model are given as below: 

£x = _xpf£ Lv f dc_ } (21) 

dt 2 8 m \ dt J 


Xp£Cd_y( ±_ w 

dt 2 8m V tit y ) 


( 2 . 2 ) 


d 2 z _ 7tpd 2 C d v 

dt 2 8 m 



(2.3) 
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where x denotes range, y denotes height, z denotes drift and W x , W y , W 2 respectively, 
the x, y and z components of the wind velocity W . This model, however, does not 
account for the effect of spin of the shell and thus fails to predict the drift due to yaw of 
repose moment. Also, it neglects the lift forces acting on the shell. 

2.3.2 Modified Point Mass Model (MPMM) 

The modified point mass model is also known as four degree-of-ffeedom model 
(three spatial degrees-of-freedom plus axial spin). Its basis is a conventional point mass 
model, in addition, the instantaneous equilibrium yaw is calculated at each time step 
along the trajectory so as to provide estimates of yaw, drag, and drift resulting from the 
yaw of repose. The equations of motion for modified point mass model are given as 
below: 



(Drag) 

(Lift) 

(Gravity and Coriolis) (2.4) 


dp _ npd 4 

~dt~ 167, 


PvCip 


a r = 


-Spl Vx[duldt ) 


(Spin damping) (2.5) 


(Yaw of repose) (2.6) 


The quantities in the above equations are defined as follows: 
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u= — , V=u-W, r=x-R 

dt 

i? =(0, — i?, 0) , [W Denotes the wind vector] 

co =(— Q cos[latitude] cos [azimuth], -Q sin [latitude], Q cos[latitude] sin[azimuth] ) 


where, Q = 7.29x1 O' 5 rad/s. (rotation of the earth), 

R = 6370320 m (radius of the earth), 

g 0 = 9.80665 [1 - 0.0026373 cos (2xlatitude) + 0.0000059 [cos (2x latitude)] 2 ], 

the axis system is as for the point mass model with x being along the line of fire. 

As may be seen from the above equations, the aerodynamic coefficient input 
required is extensive and accuracy of these aerodynamic coefficients is crucial for 
reliable estimates of range. However, for most artillery shells, the reliability of estimated 
values of aerodynamic coefficients is not high enough to inspire confidence in resulting 
range estimates from this model. In spite of this limitation, this model has the capability 
of providing reasonable accounts of end point data. It can estimate wind corrections and 
can accept non-linear aerodynamic inputs if required. 


2.3.3 Six Degree of Freedom model (SDF) 

As a final step, the most sophisticated trajectory model developed was the six- 
degree-of-freedom model having the spatial degrees of freedom, yaw degree of freedom 
in two planes and the spin. In this model there are no assumptions concerning linearised 
aerodynamics or projectile symmetry. However, the indeterminability of many of the 
initial conditions and aerodynamic coefficients which are required as input frequently 
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results in the model not giving significantly better end results than the modified point 
mass model. Thus its usage might not be justified for routine fire control work. It is 
nevertheless a powerful tool for the ammunition designer. The equations of motion for 
six-degree-of-freedom model are given as below: 


u = ( qs/m)C x - qw + rv - gs'mO + (Th/ m) (2.7) 

v= (qs/m)C v -ru + pw+ gsin</> cos 0 (2.8) 

w= ( qs/m)C : - pv + qu + g cos^cos# (2.9) 

p={qsdC,+qr(l y -l 2 )}/l , (2.10) 

q={qsdC.+rp(I,-l,)}/l y (2.11) 

f={qsdC,+pq(l,-I y )}/ 1 2 (2.12) 

<f>= P + <7 tan^sin^+r tan&cos<f> v (2.13) 

6= qcos(f>~rsm<j> (2.14) 

\jr- rcos^sec#+#sin^sec# (2.15) 


To derive the spatial position equations, above equations were used to transform 
the body-axis velocity (u, v, w) into earth fixed axis. The equations are: 

X=ucosi// cosO + v (cos y/ sin 9 sin <f>- sin \j/ cos <j>) ^ 

+ w (cos^sin#cos^+sin^sin^) +IF,. 
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(2.17) 


r=wsin^cos# + v (sin y/ sin 0 sin <f> + cos y/ cos <f>) 
+ w (sin y/ sin 0 cos <j>- cos y/ sin <j>) +W y 

Z=usw0 - v cos0sin<f>-wcos0cos0 +W Z 

Aerodynamic model used in this analysis is as given below: 
Cx=- C„ cos acosfi+C L sma- C wind cos a sin fi 

Cy=Cy» ind cos fi-Com sw.fi 

Cz= -C L cos a -Cj^M sin a cos fi-C^ sin a sin fi 

where, 

C L =C L0+ C La +C Lq (qd/2V) 

^ ywind c„+c ye p+c, r (rdl: 2V) 

C Dwind ~^d 

c, = c„ (pd/ 2V) 
c.=c. 0+ c„ a +c„( ? rf/2r) 

C„ = C M +C, f ll+C, T (rdl2V) 

Wind model used in this analysis is as given below: 

u'=u-W x cosO cos ys -W y cos0 stay/ + W z sm0 ' 

v' = v-W x (cos y/ sin 0 sin <f>- sin y/ cos <f) 

-W y (cos yr cos ^ + sin ^ sin (9 sin 
-R^cos^sin^ 


(2.18) 
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w'=w-W x (cos y/ sin 9 cos sin y/ sin </>) 

— W (sin y/ sin 6 cos <f>- cos y/ sin (f) 
-W z cos <9 cos ^ 


thus, a=tan 1 -^7 and /?=tan 1 

\u ) \u J 

where W x , W y , W z are the wind components blowing toward the 
x, y, and z directions. 


2.4 Generation of simulated flight data 

The capability of the proposed schemes to estimate parameters from spatial 
coordinates of the artillery shell is investigated using the following sets of flight data 
generated through the above models: 


SI: In this case the flight data (t,x,y) were generated through PMM for 
different elevations using C d varying with Mach number as listed in Table 2.1. Flight 
data thus obtained were used to estimate equivalent value of C d and its dependence on 
Mach number respectively. 

S2: In this case the flight data (t,x,y) were generated using C d varying with 

Mach number for elevations 10, 20 so on up to 60 Deg. respectively and only end points 
(where the projectile strikes the ground) had been recorded for estimation purpose. In 
most of the field trails, radar tracked data may not be easily available. Conventionally 
only the end point data are measured and recorded. Hence to investigate the suitability of 
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the proposed schemes, these were applied on end point data to estimate an equivalent 


value of C d . 

S3: In order to check robustness of the proposed schemes while dealing noisy 
data, flight data were corrupted with known noise. Thus, in this case simulated pseudo 
random noise of varying intensity was added to the flight data that had been generated in 
SI which intern was used to estimate the functional dependence of Mach number on C d . 

A measure of confidence level was estimated through the estimation of Cramer-Rao 
bound. 


After getting some initial success, the simulation model to generate flight data 
was progressively made exhaustive by incorporating more number of aerodynamic 
parameters. 3-Dimensional flight data of an artillery shell were generated using MPMM 
and SDF models. 


S4: In this case the flight data (t,x,y ,z) were generated through MPMM using 
constant value of parameters C d , C La , C ma and C lp . In this model number of parameters 

to be estimated are four. These data were used to check suitability of proposed scheme in 
handling more then one parameter for estimation. 

S5: In this case the flight data (t,x,y,z) were generated through MPMM by 

modeling the variation of aerodynamic parameters (C d ’ C La » C m a and with Mach 

number. Flight data thus obtained were used to estimate equivalent values of parameters 
and Mach number dependence of all the parameters. 
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S6: In this case the flight data (t, x,y ,z) were generated through SDF model 


using constant value of parameters C d , C La , C ma and C lp . In this model number of 

parameters to be estimated are still more. So, in order to check the applicability of the 
proposed schemes, these data were used to test whether the proposed schemes can 
estimate all these parameters at one time. 

S7 : In this case the flight data (t,x,y ,z) were generated through SDF model by 

modeling the variation of aerodynamic parameters (C d , C Ia , C ma and C, p ) with Mach 

number. Flight data thus obtained were used to estimate equivalent values of parameters 
and Mach number dependence of all the parameters. 
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CHAPTER 3 

RESULTS AND DISCUSSIONS 


3.1 General 

In this chapter results of the estimated parameters will be presented for artillery 
shell utilizing the schemes proposed for flight data in chapter 2. Before presenting the 
detailed results of estimated parameters and discussing their accuracy, a brief 
discussion is given below about the initial values, number of data points, sampling rate 
and intensity of measurement noise. 

Different combinations of initial guess values of the parameters were used to 
start the iterative algorithm of Maximum Likelihood (ML) method. In general, results 
were found not much affected by the choice of the initial values except for a few cases. 
It is therefore decided to present results for typical sets of initial values respectively. 

The optimal number of data points chosen varies from case to case as will be 
explained accordingly. The sampling rate had been fixed at t - 0.01 sec for almost all 
the cases. The simulated flight data, as explained in chapter 2, in terms of spatial 
coordinates (x ,y,z) had been generated and used for estimation. 

A study was also carried out to see the effect of measurement noise on the 
estimation of parameters. Accordingly, simulated pseudo random noise of varying 
intensity was added to the flight data. The noise was simulated by generating 
successive uncorrelated pseudo random numbers having normal distribution with zero 
mean and assigned standard deviation. The intensity of the noise differs from case to 
case and results are presented for the same. Cramer-Rao bounds had been calculated for 
both no noise and noise cases and they seem to be well within the limits. This shows 
high confidence level of the estimated parameters. 
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3.2 Estimation of equivalent parameters using scheme 1 

In this scheme motivation, as discussed earlier, was to estimate equivalent 
parameters which can be assumed to be an approximation for the variation of 
parameters with Mach number. This scheme was applied to various simulated flight 
data generated in order to estimate equivalent parameters as discussed below: 

3.2.1 Estimation of equivalent C d 

In order to test the code, flight data were generated using PMM with a constant 
value of C d . Using this simulated flight data the value of aerodynamic parameter Q 
was estimated using ML method. The estimated value of C d exactly matched with that 
of the true value. Table 3.1 presents this result. 


Table 3.1 Comparison of estimated and true values of C d ; scheme 1. 


SL. No. 

Parameter 

True Value 

Initial Guess 
Value 

Estimated 

Value 

1 . 

c d 

0.30 

0.50 

Or 

0.10 

0.30 

(0)’ 


*Cramer-Rao bound. 


After validating the code by above approach, equivalent value of C d was 
estimated using the simulated flight data SI. For this purpose, flight data were 
generated for a set of different elevations and an average value of C d , pertaining to 
each elevation, was estimated using scheme 1 as shown in Table 3.2 respectively. 

sfTo A • •** ■" X&hx- - 
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Table 3.2 Equivalent value of C d estimated; scheme 1, flight data SI. 


SL. No. 

Elevation (0), degrees 

Equivalent C d 

1 . 

10 

0.267803 

(0.0000026)’ 

2. 

15 

0.274797 

(0.0000015) 

3. 

20 

0.276251 

(0.0000010) 

4. 

25 

0.275646 

(0.00000077) 

5. 

30 

0.274519 

(0.00000061) 

6. 

35 

0.273226 

(0.00000050) 

7. 

40 

0.271675 

(0.00000042) 

8. 

45 

0.270041 

(0.00000037) 

9. 

50 

0.268703 

(0.00000033) 

10. 

55 

0.267809 

(0.0000003) 

11. 

60 

0.267432 

(0.00000028) 

12. 

65 

0.267479 

(0.00000026) 

13. 

70 

0.267682 

(0.00000026) 

*Cramer-Rao 1 

sound 


Average equivalent value of C d i.e. 0.271005 was obtained by taking average 
of individual estimated C d values. Using this equivalent C d value, estimated flight data 
(t,x,y) were compared with that of measured data for different elevations as shown in 
Fig. 3.1. The error in both range and height between measured and estimated responses 
are plotted with range as shown in Fig. 3.2. Table 3.3 compares the range between 
measured and estimated responses for different elevations. Comparing measured and 
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Fig 3.1 Comparison of measured and estimated trajectories obtained using 
equivalent C d ; scheme 1, flight data SI. 
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Fig 3.2 Comparison of error in range and height between measured and estimated 
trajectories obtained using equivalent C d ; scheme 1, flight data SI. 
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estimated range, it can be observed that the maximum errors in range and height are 
well within one percent of the maximum range/height . Further low values of Cramer- 
Rao bound show higher level of confidence in estimation. 


Table 3.3 Comparison of measured and estimated range using equivalent C d ; 
Scheme 1, flight data SI. 


SL. No. 

Elevation (0), 

degrees 

Measured Range, m 

Estimated Range, m 

1 . 

10 

12193.078 

12131.827 

2. 

32 

i 

21100.430 

21229.985 

3. 

52 

23088.499 

22773.323 

4. 

70 

16701.777 

16403.438 


3.2.2 Estimation of equivalent C d using only end point data 

In this case, the flight data were generated as stated in S2. Scheme 1 was 
applied to estimate equivalent value of C d using this data. The estimated value of 
was found to be 0.269074. Using this equivalent C d value, estimated flight data 
(t, x,y) were compared with that of measured data for different elevations as shown in 
Fig. 3.3. The error in both range and height between measured and estimated responses 
is also plotted with range as shown in Fig. 3.4. Table 3.4 compares the range between 
measured and estimated responses for different elevations. Notwithstanding the 
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Fig 33 Comparison of measured and estimated trajectories using 
equivalent C d ; scheme 1, flight data S2. 
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Fig 3.4 Comparison of error in range and height between measured and estimated 
trajectories using equivalent ; scheme 1, flight data S2. 
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difference in range, this scheme could be very useful in getting some average value of 
C d for a particular type of shell through end point data. 


Table 3.4 Comparison of measured and estimated range using equivalent C d ; 
scheme 1, flight data S2. 


SL. No. 

Elevation (<9), 

degrees 

Measured Range, m 

Estimated Range, m 

1 . 

10 

12193.078 

12170.877 

2. 

25 

19038.396 

19329.929 

3. 

45 

23219.449 

23169.223 

4. 

70 

16701.777 

16493.113 


3.2.3 Estimation of equivalent values of C d ,C La , C ma and C lp from S5 

In order to test the code flight data were generated as stated in S4. ML method 
was applied to estimate these parameters. It was found that the estimation algorithm 
failed to estimate C la and C ma . However other two parameters were well estimated. 

Close look at the Eqn. (2.4) reveals that C La and C ma were highly correlated as far as 
estimation is concerned. Hence, it was decided to fix one of these parameters. Thus 
C Ia was Jcept fixed to its tpje value during estimation. Table 3.5 presents the estimated 

results, which are exactly same as true values used for generation of flight data. 
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Table 3.5 Comparison of estimated and true value of parameters; 
scheme 1, flight data S4. 



Parameter 

True 

Value 

Initial Guess 

Value 

Estimated 

Value 

1 . 


0.30 

0.5 

0.30 

(0)’ 

2. 

Cu, 

(Fixed) 

2.0 

2.0 

2.0 

3. 

c 

mcc 

3.0 

4.0 

3.0 

(0) 

4. 

C . P 

-0.03 

-0.05 

-0.03 

(0) 


*Cramer-Rao bound 


After validating the code by above approach, equivalent values of the 
parameters were estimated using scheme on S5. For this purpose, flight data were 
generated for a set of different elevations and average values of the parameters, 
pertaining to each elevation, were estimated. In the estimation process was kept 
fixed to 2.0. Table 3.6 lists the average values of the parameters estimated for different 
elevations. Table 3.7 compares measured and estimated range, drift. Referring Table 
3.6 it can be observed that roll damping parameter (C lp ) exhibit large scatter. This is 

expected as C lp for shell is a weak derivative (parameter) and its scatter does not effect 

strongly on range or drift. Notwithstanding this scatter, the comparison of estimated 
and measured flight data shows excellent matching for all practical purposes. 
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Table 3.6 Equivalent values of parameters estimated; scheme 1, flight data S5. 


SL. No. 

Parameter 

o 

o 

II 

o 

O 

eo 

II 

0 = 45° 

1. 

Q 

0.26810 

(0.0000026) 

0.27415 

(0.0000006) 

0.26954 

(0.0000003)* 

2. 

c* 

(Fixed) 

2.0 

2.0 

2.0 

3. 


1.90432 

(0.0068) 

1.45614 

(0.00094) 

1.49519 

(0.0006) 

4. 


-0.07802 

(0.00086) 

-0.15129 

(0.0001) 

-0.15915 

(0.00006) 


*Cramer-Rao bound. 


Table 3.7 Comparison of measured and estimated range, drift using equivalent 
parameters; scheme 1, flight data S5. 


SL. 

Elevation, 

Range, m 

Drift, m 

No. 

e 

Measured 

Estimated 

Measured 

Estimated 

1 . 

10° 

12088.44 

12092.48 

89.51 

91.12 

2. 

30° 

20390.55 

20376.29 

386.46 

383.873 

3. 

45° 

23002.13 

22866.49 

690.06 

677.96 


3.2.4 Estimation of equivalent values of C d , C La , C ma , C yfi , C nfi andC, p from S7 

In order to test the code flight data were generated as stated in S6. ML method 
was applied to estimate these parameters. It had been observed that all the four 




























parameters were estimated together and there values were exactly equal to their true 
values as shown in Table 3.8. There had been no correlation problem observed here. 
Due to symmetry of the shell following modeling was incorporated in estimation 
algorithm i.e. C yP = -C La and C nf =-C ma . 

Strangely enough, the estimation of C mq posed a serious problem. During 
estimation iteration, the value of C mq was changing sign abruptly from -ve to +ve. At 
this point the trajectory computation failed to yield any reasonable result. For shell it 
was found that C mq did not effect range or drift appreciably. Hence it was decided to 

omit C mq from estimation algorithm. 

Table 3.8 Comparison of estimated and true value of parameters; 
scheme 1, flight data S6. 


SL.No. 

Parameter 

True 

Value 

Initial Guess 

Value 

Estimated 

Value 

1. 


0.30 

0.5 

0.30 

(0)* 

2. 

C. 

2.0 

3.0 

2.0 

(0) 

3. 


3.0 

4.0 

3.0 

(0) 

4. 

C ip 

-0.03 

-0.05 

-0.03 

(0) 


*Cramer-Rao bound 


After validating the code by above approach, equivalent values of the 
parameters were estimated using scheme on S7. However, this scheme foiled to 
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estimate equivalent values of the parameters to any reasonable accuracy when full 
trajectory data were used. When small part of the flight data (initial 2 sec) was used 
equivalent values of all the parameters were reasonably estimated. Further analysis is 
required to investigate the failure of estimation when applied for full trajectory. 

3.2.5 Estimation of constant C d from noisy flight data 

In order to investigate the robustness of the present scheme in handling noisy 
data, SI was generated with constant C d and was corrupted with known noise. Using 
this noisy data once again constant C d was estimated. Table 3.9 shows the comparison 
of constant C d estimated for different noise levels. Observing the values of 
C d estimated along with Cramer-Rao bound, it could be concluded that proposed 
scheme is robust enough to process noisy data to an acceptable level of confidence. 

Table 3.9 Comparison of estimated and true values of C d ; scheme 1. 


SL.No. 

True Value 

C d Estimated with different noise levels 

5 m 

10 m 

50 m 

100 m 

1 . 

0.30 

0.3000004 

(0.000002)* 

0.3000008 

(0.000005) 

0.3000043 

(0.000026) 

0.3000087 

(0.000053) 

*Cramer- 

Rao bound 


3.3 Estimation of Mach number dependence of aerodynamic 
parameters using scheme 2 

In this scheme in order to capture the Mach number dependence of aerodynamic 
parameters the modeling in estimation algorithm was done using series approximation 
of the parameters as a function of Mach number and the constants of this series were 
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estimated applying ML method. 

3.3.1 Estimation of C d vs Mach number from SI. 

In this case flight data SI were used to estimate C d varying with Mach number. 
In this scheme, in estimation algorithm C d was expressed as a polynomial function of 
Mach number as C d = aO + al * M + al * M 2 . The constants aO , al and al have been 
estimated through ML method. In order to capture the whole range of Mach number, 
flight data were generated at different initial velocities of 818m/s and 360m/s at an 
elevation of 7 Deg. . From the flight data it was observed that the shell when fired at 
818m/s covers Mach number range from 2.404 to 1.159. Further when fired at 360m/s, 
the shell covered Mach number from 1.058 to 0.904. Thus entire range of Mach number 
was covered. 

The values of a0,al and al estimated for both the charges are presented in 
Table 3.10. The values of estimated C d Vs M is given in Table 3.11 for both the 

charges and is compared with true C d Vs M in Fig 3.5. 


Table 3.10 Estimated values of constants; scheme 2, flight data SI. 


SL. No. 

V = 818 m/s 

V =360 m/s 

Parameter 

Value 

Parameter 

Value 

1 . 

aO 

-0.563540 

(0.00229) 

aO 

-3.326722 

(0.67713)* 

2. 

al 

-0.235335 

(0.00257) 

al 

5.320928 

(1.38287) 

3. 

i 

al 

0.041198 

(0.00069) 

al 

-1.694348 

(0.70489) 


*Cramer-Rao bound 






















Table 3.11 Estimated C d Vs M; scheme 2, flight data SI. 


SL. No. 

Charge (V) = 818 m/s 

Charge (V) =360 m/s 

M 

c d 

M 

Cd 

1. 

2.359184 

0.237639 

1.035000 

0.365411 

2. 

2.309184 

0.239790 

1.030000 

0.356300 

3. 

2.259184 

0.242146 

1.025000 

0.347105 

4. 

2.209184 

0.244708 

1.020000 

0.337825 

5. 

2.159184 

0.247477 

1.015000 

0.328460 

6. 

2.109184 

0.250451 

1.010000 

0.319011 

7. 

2.059184 

0.253631 

1.005000 

0.309477 

8. 

2.009184 

0.257018 

1.000000 

0.299858 

9. 

1.959184 

0.260610 

0.995000 

0.290154 


1.909184 

0.264408 

0.990000 

0.280366 

11. 

1.859184 

0.268413 

0.985000 

0.270493 

12. 

1.809184 

0.272623 

0.980000 

0.260536 

13. 

1.759184 

0.277039 

0.975000 

0.250493 

14. 

1.709184 

0.281661 

0.970000 

0.240366 

15. 

1.659184 

0.286490 

0.965000 

0.230154 1 

16. 

1.609184 

0.291524 

0.960000 

0.219858 

17. 

1.559184 

0.296764 

0.955000 

0.209477 

18. 

1.509184 

0.302210 

0.950000 

0.199011 

19. 

1.459184 

0.307862 

0.945000 

0.188460 

20. 

1.409184 

0.313721 

0.940000 

0.177824 

21. 

1.359184 

0.319785 

0.935000 

0.167104 

22. 

1.309184 

0.326055 

0.930000 

0.156299 

23. 

1.259184 

0.332531 

0.925000 

0.145410 

24. 

1.209184 

0.339213 

0.920000 

0.134436 

25. 

1.159184 

0.346101 




3.3.2 Estimation of C d vs Mach number from S7 

In order to study the effect of sensitivity of trajectory modeling in estimation 
algorithm, the estimation algorithm only had PMM, whereas the flight data were 
generated was using complete SDF. It was observed that as far as C d vs M variation is 
concerned, s imp le PMM is adequate enough to capture this variation. 
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Fig 3.5 Comparison of C d Vs M between estimated and measured values; 
scheme 2, flight data SI. 



Fig 3.6 Comparison of C d Vs M between estimated and measured values; 
scheme 2, flight data S7. 
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The values of aO,al and al estimated for both the charges are given in 
Table3.12. The values of estimated C d Vs M is given in Table 3.13 for both the charges 
and is plotted with true C d Vs M in Fig 3.6. 


Table 3.12 Estimated values of constants; scheme 2, flight data S7. 


SL. No. 

V = 818 m/s 

V =360 m/s 

Parameter 

Value 

Parame 

ter 

Value 

1 . 

a 0 

-0.548640 

(0.00228) 

a0 

-5.449019 

(0.6611)* 

2. 

al 

-0.212708 

(0.00256) 

al 

9.736826 

(1.3507) 

3- 

al 

0.034592 

(0.000690) 

al 

-3.982691 

(0.6887) 


*Cramer-Rao bound 


3.3.3 Estimation of C d , C La , C ma and C lp vs Mach number 

Application of scheme 2 to capture Mach number dependence of ^La ’ C ma 5 C !p 
together with C d posed numerical difficulties. It is recommended that these parameters 
needed to be expressed using appropriate series expansion. Simple series 
approximation as was used in case of estimating C d vs M was not suitable here. 

3.4 Estimation of Mach number dependence of aerodynamic 
parameters using scheme 3 

In this scheme in order to capture the Mach number dependence of aerodynamic 
parameters the estimation was done by splitting the whole trajectory into different sets 
of points and average values of the parameters pertaining to each set of points were 
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estimated. The estimated parameters correspond to average Mach number of the set of 
points considered for estimation. 


Table 3.13 Estimated C d Vs M; scheme 2, flight data S7. 


SL. No. 

Charge (V) = 818 m/s 

Charge (V) =360 m/s 

M 


M 

c d 

1 . 

2.354109 

0.239605 

1.037632 

0.36614 

2. 

2.304109 

0.242184 

1.027632 

0.351025 

3. 

2.254109 

0.244935 

1.017632 

0.335113 

4. 

2.204109 

0.24786 

1.007632 

0.318405 

5. 

2.154109 

0.250957 

0.997632 

0.3009 

6. 

2.104109 

0.254227 

0.987632 

0.282599 

7. 

2.054109 

0.257671 

0.977632 

0.263501 

8. 1 

2.004109 

0.261287 

0.967632 

0.243606 

9. 

, 1.954109 

0.265076 

0.957632 

0.222915 

10. 

1.904109 

0.269039 

0.947632 

0.201428 

11. 

1.854109 

0.273174 

0.937632 

0.179144 

12. 

1.804109 

0.277482 

0.927632 

0.156063 

13. 

1.754109 

0.281963 

0.917632 

0.132186 

14. 

1.704109 

0.286617 



15. 

1.654109 

0.291444 



16. 

1.604109 

0.296444 



17. 

1.554109 

0.301617 



18. 

1.504109 

0.306963 



19. 

1.454109 

0.312482 



20. 

1.404109 

0.318174 



21. 

1.354109 

0.324038 



22. 

1.304109 

0.330076 



23. 

1.254109 

0.336287 




24. 

1.204109 

0.342671 



25. 

1.154109 

0.349227 




3.4.1 Estimation of C d \ s Mach number 

In this case the flight data SI was used to estimate C d varying with Mach 
number. In this scheme the whole trajectory (flight data) was split into different sets of 
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100 points each. Then an average value ofQ, pertaining to each set of points, was 
estimated by processing these sets of points. The initial conditions (velocity and 
elevation (0)) required for generation of estimated response for each set of points was 
calculated by differentiating x and y with respect to time. The estimated C d Vs M is 
plotted against original C d Vs M as shown in Fig 3.7 for different elevations and the 

respective values, for an elevation of 10 Deg., are shown in Table 3.14. It can be seen 
from Table 3.14 and figures that the estimated values matched closely with that of the 
true values. 

3.4.2 Estimation of C d ,C La ,C ma and C lp vs Mach number 

In order to estimate the functional dependence of all the parameters with Mach 
number we require all the initial conditions of the motion variables. As the estimation 
algorithm was made exhaustive, using MPMM and SDF models, the number of initial 
conditions required became large. It was not possible to get all these initial conditions 
through spatial coordinates (x,y,z) only. This was only possible when estimation 

algorithm used PMM. 

3.4.3 Estimation Of Q Vs Mach number from noisy data 

In this case the flight data generated in S3 (noisy data) was used to estimate C d 
varying with Mach number. To apply this scheme, for estimating QVs M, 
differentiation of x and y were essential. For noisy data differentiation could not be 
done to yield reasonable result as presented in Fig 3.8. To overcome this difficulty, the 
noisy data was smoothened by fitting a polynomial. 
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Mach Number ( M ) 


Mach Number ( M ) 



Fig. 3.7 Comparison of C d vs M between estimated and measured values; 
scheme 3, flight data SI. 








Table 3.14 Estimated C d Vs M; scheme 3, flight data S7. 


SL. No. 

Mach Number 

c d 

1 . 

2.340629 

0.237474 

(0.002630)* 

2. 

2.226243 

0.242646 

(0.002950) 

3. 

2.121262 

0.248222 

(0.003291) 

4. 

1.934627 

0.260397 

(0.004040) 

5. 

1.851056 

0.267130 

(0.004449) 

6. 

1.773072 

0.273649 

(0.004882) 

7. 

1.700093 

0.280969 

(0.005339) 

8. 

1.631685 

0.287871 

(0.005820) 

9. 

1.567552 

0.294438 

(0.006323) 

10. 

1.450988 

0.306815 

(0.007394) 

11. 

1.398073 

0.312466 

(0.007959) 

12. 

1.301720 

0.324071 

(0.009138) 

13. 

1.257749 

0.331029 

(0.009751) 

14. 

1.216238 

0.338202 

(0.010378) 

15. 

1.177040 

0.345416 

(0.011018) 

16. 

1.104875 

0.362316 

(0.012331) 

17. 

1.071931 

0.366093 

(0.012994) 

18. 

1.041634 

0.364824 

(0.013639) 

19. 

1.014675 

0.349604 

(0.014232) 

20. 

0.992291 

0.317082 

(0.014728) 


* Cramer-Rao bound. 
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Fig 3.8 Comparison of C d Vs M between estimated and measured values; 
scheme 3, flight data S3. 



Fig 3.9 Comparison of C d Vs M between estimated and measured values 
after smoothening; scheme 3, flight data S3. 
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For an elevation of 10 Deg. the noisy data is approximated using a 6 th order 
polynomial given by: 

Range: -3.6e-07*x 6 +0.000 13* x 5 -0.011** 4 + 0.513** 3 

-20.05** 2 +804.12** + 0.718. ' 

Height: - 1.47e-06*/ + 0.00013*/ -0.0058*/ + 0.193*/ 

, ( 3 - 2 ) 

-8.57*/ + 142.15*7-0.103. 

The flight data, smoothened using the polynomial approximation, was used for 
estimation. The estimated C d Vs Mach number compares well with that of true QVs 

M as shown in Fig 3.9. For future work, flight data may be processed using standard 
filtration techniques before being used for parameter estimation. 

3.5 Estimation of Mach number dependence of aerodynamic 
parameters using scheme 4 

3.5.1 Estimation of C d , C la , C ma and C lp vs Mach number 

In order to estimate the dependence of all the parameters at a time with Mach 
number, the earlier schemes were inadequate. In order to avoid such difficulties an 
alternate scheme have been proposed for field application. 

Flight data were generated corresponding to different charges (velocities) and 
the initial part of the trajectory was used for estimation. In field application this 
amounts to firing few rounds of different charges at some elevation. Initial 200 points 
of the trajectory were taken to estimate average values of the parameters (C d , C La , 
C , C, ) for different charges respectively. Hence the equivalent parameters 
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estimated belong to the average Mach number for the set of 200 points. Thus for 
different charges fired we get equivalent parameters corresponding to that Mach 
number. The aerodynamic parameters thus estimated as function of Mach number are 
presented in Fig 3.10. Referring Fig 3.10, it can be observed that for most of the 
parameters, the Mach number dependence was well estimated. However estimated 
value of C lp showed large scatter. 
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Mach Number (M) Mach Number (M) 


Fig 3.10 Comparison between measured and estimated parameters; 
scheme 4, flight data S7. 












CHAPTER 4 


CONCLUSIONS AND FUTURE SUGGESTIONS 

4.1 Conclusions 

In the present work, schemes were proposed to estimate, both equivalent and 
Mach number dependence of aerodynamic parameters from radar tracked flight data 
(x,y,z) of an artillery shell. Four schemes using ML method are proposed for parameter 
estimation. The proposed schemes show promise to estimate the aerodynamic parameters. 

Scheme 1 is used to estimate equivalent aerodynamic parameters and scheme 2 to 
4 are used to estimate functional dependence of parameters on Mach number. Scheme 2 
assumes series approximation of the parameters as a function of Mach number in 
estimation algorithm. In scheme 3, whole trajectory was split into different sets of points 
and average values corresponding to each set of points were estimated. In scheme 4 the 
flight data was generated at a particular elevation corresponding to different charges and 
small initial part of the trajectory was used to estimate the estimate the Mach number 
dependence of parameters. 

Based on all the schemes it can be concluded that scheme 1 can be suitably 
applied to get average values of aerodynamic parameters. The differences in range/drift 
using the average values are well within 1% of the maximum range/drift. Scheme 2 and 3 
are best suited for capturing the Mach number dependence of C d . However scheme 3 
faced difficultie s in order to estimate the Mach number dependence of all the parameters, 
scheme 4 was best suited for this. 
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4.2 Suggestions for future work 


1. In the present work, simulated flight data were used to estimate aerodynamic 
parameters of the artillery shell due to non-availability of the real flight data. In order to 
enhance the capability of the proposed schemes, these should be validated using real 
flight data. 

2. To deal with noisy data better filtration techniques should be used to pre process the 
data before used for parameter estimation. 

3. Scope of the scheme 2 may be expanded to capture Mach number dependence of all 
the par ame ters by assuming appropriate model, representing the variation of parameters 
with Mach number, in estimation algorithm. 

4. All the schemes were applied to the simulated data generated for standard atmospheric 
conditions. Applicability of the proposed schemes to estimate parameters using flight 
data under non-standard atmospheric conditions may be explored. 
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APPENDIX A 


MAXIMUM LIKELIHOOD (ML) METHOD 

The parameter estimation algorithm used in this study is the conventional 
‘Maximum Likelihood Method’. The programs were written in ‘C’ using the MLE 
algorithm. The algorithm is given as follows: 

Let a mathematical model of a system be represented by the following equations: 

X{t,) = AX(t) + BU(t,) (A.l) 

Z(t,) = CX(t) + Gn, (A.2) 

where X is the state vector, U the control vector, Z the measured (observation) 
response and «, is the noise vector respectively. The matrices A and B contain the 

unknown stability and control derivatives (parameters). If there is no state noise and the 
matrix GG is known, then the Maximum Likelihood estimator minimizes the cost 
function, 

J(#H 1/2) t {[ z (t)-Z,(t) ] r (<3G r )~' [ Z(t,)-Z e (t,) ]} (A.3) 


where GG T is the measurement noise covariance matrix and Z g (t, ) is the computed 
response estimate of Z at f, for a given value of the unknown parameter vector 6 . The 
cost function is a function of the difference between the measured and computed time 
histories. In our case, we have taken GG T = 1 for the cost function. To minimize the cost 
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function «/(#), we can apply the Newton-Raphson algorithm, which chooses successive 
estimates of the vector 9 of unknown coefficients. If k is the iteration number, then 


<V, = e„ - [ (e , ) Y [ v j j(e t ) ] , (a.4) 

where the first gradient is defined by 

V » j («) = -2{[ Z (':)- Z »('-)] r ( GG T'[ V «' Z »('.)]}. (A- 5 ) 

and the second gradient is approximated by 

= Z{[^ Z »('*)] , (GG r F[V,Z,((,)]}. (A.6) 


The matrix GG T is approximated by a diagonal matrix, with the diagonal 
elements given by. 


{[Z(>,)- Z »(t,)][ Z (>,)- Z < (<,)] r } 


s 

/=] 


iV 


(A. 7) 


where N is total number of data points. The term V s Z e (r f ) occurring in the first 

and the second gradients is evaluated by finite difference approximation for each of the 
sensitivity coefficients present in it. For example, to evaluate sensitivity coefficient 

da 

dC, 


■ , the equations of motion are solved to obtain a for C mq + A C mq and 


'mq 


mq 


mq A G mq 
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Let the respective values of a be a + and a~, then we have 


da _ a + — a' 
d C mq 2 A C mq 


The reliability of the estimated parameter in terms of the Cramer-Rao bounds is 
given by the square root of the diagonal elements of [v^ j(6 k ) ]~\ which is anyway 
evaluated as part of the ML algorithm. 

In the present work, for estimating the aerodynamic parameters of the artillery 
shell, the vector of unknown parameters is given by d~ [ C d ,C la ,C ma ,C lp ] for a 

typical case with the Z (r, ) and Z e (t , ) responses being computed by the spatial 
coordinates x,y and z respectively. 
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